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Executive Summary 

Current forecasts of growth in air traffic demand indicate that the existing capacity 
at many of the busiest airports will be inadequate to avoid exponentially increas- 
ing delays. NASA is investigating new technologies and operational tools that will 
reduce some of the current constraints on airport capacities. 

To meet its objective of assisting U.S. industry with the technological challenges 
of the future, NASA must identify research areas that have the greatest potential 
for improving the operation of the air transportation system. Therefore, NASA 
seeks to develop the ability to evaluate the potential impact of various advanced 
technologies. By thoroughly understanding the economic impact of advanced 
aviation technologies, and by evaluating how these new technologies would be 
used within the integrated aviation system, NASA aims to balance its aeronautical 
research program and help speed the introduction of high-leverage technologies. 
To accomplish this goal, NASA is building an Aviation System Analysis Capa- 
bility (ASAC). 

The ASAC is envisioned primarily as a process for understanding and evaluating 
the impact of advanced aviation technologies on the U.S. economy. ASAC con- 
sists of a diverse collection of models, databases, analysts, and other individuals 
from the public and private sectors brought together to work on issues of common 
interest to organizations within the aviation community. ASAC also will be a re- 
source available to those same organizations to perform analyses; provide infor- 
mation; and assist scientists, engineers, analysts, and program managers in their 
daily work. 

The ASAC Airport Capacity Model and the ASAC Airport Delay Model support 
analyses of technologies addressing airport capacity. There are two primary con- 
straints on airport capacity that NASA technologies address. The first is arrival 
runway occupancy time (AROT), the time from when an arriving aircraft crosses 
the end of the runway until it turns off the runway. For safety reasons, only one 
aircraft may occupy a runway at any time. Reducing AROT may allow a quicker 
tempo on the runway, thus increasing the airport’s capacity. Reducing the uncer- 
tainty in AROT may allow air traffic controllers to shorten the spacing between 
successive arrivals without compromising safety, also increasing airport capacity. 


iii 



The amount of time an aircraft requires on a runway and its standard deviation 
(uncertainty) are input values to the ASAC Airport Capacity Model. The user can 
compare airport capacity with typical AROT values (which are airport specific) 
versus capacity with reduced AROT values anticipated as a result of a new tech- 
nology. 

The second major factor affecting airport capacity is the separation that must be 
maintained between any two airborne craft in the terminal environment. The cur- 
rent separation standards are determined on the basis of safety concerns. Some 
NASA initiatives provide technology that may enable aircraft to operate in closer 
proximity without compromising safety. 

Separation requirements are also ASAC Airport Capacity Model inputs. The user 
may set the separation criteria (which are weather dependent) to those appropriate 
to a new technology and examine the resulting capacity impact. 

Other Federal Aviation Administration operational rules may become obsolete if 
new technologies are developed. The user can direct the ASAC Airport Capacity 
Model to either take into account or disregard these rules, as appropriate to the 
technology under examination. 

The ASAC Airport Delay Model allows the user to estimate the minutes of arrival 
delay for an airport, given its (weather dependent) capacity. Historical weather 
observations and demand patterns are provided by ASAC as inputs to the delay 
model. The ASAC economic models can translate a reduction in delay minutes 
into benefit dollars. 

The ASAC models have been developed to account for the peculiarities of each 
modeled airport, determined on the basis of interviews with air traffic controllers 
responsible for operations at the airport. The results of each airport model, using 
current AROTs and separations, have been reviewed by these air traffic control- 
lers, and their concurrence with the models’ results has been obtained. 
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Chapter 1 

The ASAC Airport Capacity Model 


NASA’s Aviation System Analysis Capability (ASAC) Airport Capacity Model 
estimates the capacity of an airport as a function of weather. Federal Aviation 
Administration (FAA) procedures, traffic characteristics, and the level of technol- 
ogy available. Airport capacity is presented as a Pareto frontier of arrivals per 
hour versus departures per hour. This frontier is the set of points at which arrival 
rate and departure rate cannot be simultaneously increased. 

Figure 1-1 gives an example capacity curve. The figure indicates that, when de- 
partures are given priority, Newark International Airport can accommodate up to 
57 departures per hour. The airport can integrate up to 46 arrivals per hour into the 
departure stream while maintaining that departure rate. The airport can increase 
the arrival rate to more than 46 per hour only by decreasing the departure rate. By 
decreasing the departure rate, the airport can increase the maximum arrival rate to 
up to 50 per hour. (An airport capacity curve is not necessarily made up of 
straight-line segments as in the example.) 

Figure 1-1. Example of Airport Capacity: 

Newark International Airport 



Actual airport capacity varies with, among other factors, ceiling, visibility, wind 
speed, and wind direction, so that a complete specification of airport capacity is a 
family of curves like that in Figure 1-1. The ASAC Airport Capacity Model pres- 
ents a capacity curve for each of five weather conditions: visual approaches 
(VMC1); radar approaches with visual completion (VMC2) and three categories 
of instrument approaches — IMC1, IMC2, and IMC3. The airport capacities may 
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vary in the different IMC cases, as some runways may not be available as the 
weather deteriorates. The ceiling and visibility criteria that define weather condi- 
tions and runway availability were obtained from air traffic controllers at each 
modeled airport and from published information about these airports. 

The capacity models present the weather-dependent set of frontiers for each air- 
port configuration. An airport configuration is the specification of geometry and 
strategy; that is, which runways are being used and how they are being used. For 
example, a typical configuration at Detroit Metropolitan Wayne County Airport 
(DTW) is to use runways 21L, 21C, and 21R (geometry), with 21L used for arri- 
vals only, 2 1C for departures only, and 21R used for departures and arrivals in the 
proportion needed to best satisfy demand (strategy). 

Although each airport model has unique configuration logic, all the airports use 
the same analytical tools to translate technology, traffic characteristics, and opera- 
tional rules into runway capacity. This suite of tools comprises the Runway Ca- 
pacity Module, discussed in the next section. 

The Runway Capacity Module 

The Runway Capacity Module analyzes runway capacity, incorporating parame- 
ters related to the separation requirements for aircraft, runway occupancy times, 
controller information, and operational rules that advanced technologies can af- 
fect. The Runway Capacity Module takes a controller-based view of airport op- 
erations. The limitations on the quality of information accessible to the 
controller — such as aircraft position and speed — directly affect the spacing re- 
quired for safe operation of aircraft streams. 

Table 1-1 lists the key parameters used in the AS AC Runway Capacity Module. 
This module determines controller-applied separation of aircraft at the beginning 
of the common approach path (the last point at which instructions may be issued 
to the aircraft) on the basis of aircraft performance parameters, required separa- 
tion, and the requirement that two arriving aircraft not occupy the runway at the 
same time. Uncertainties about aircraft position and speed require the controller to 
apply separations greater than the minimum required at the beginning of the 
common path, to achieve the required separation from the outer marker to the 
runway threshold with a specified level of statistical confidence. Uncertainty 
about runway occupancy time similarly requires the controller to apply greater 
separation to ensure that two aircraft do not occupy the runway at the same time. 

To model a particular technology, the user may toggle (on or off) or set parame- 
ters in specific operational rules affecting capacity, as appropriate. One such rule 
is the required departure hold time after a “heavy” jet departure. Current FAA 
procedures classify an aircraft as heavy if its gross takeoff weight can exceed 
255,000 pounds. Such aircraft generate strong wingtip vortexes. The current rule 
requires that, after a heavy jet departure, the next departure on the same or parallel 
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The ASAC Airport Capacity Model 


runway within 2,500 feet cannot occur for 2 minutes. A related rule specifies a 
3-minute wait before allowing an arriving or departing aircraft to fly through the 
flight path of a heavy jet at a runway intersection. The ASAC user may leave the 
wake vortex rules as is, adjust either or both hold times, or toggle the rules off. 

Table 1-1. Capacity Model Parameters 


Symbol 

Definition 

Pi 

Fraction of aircraft in class / 

S,j 

Miles-in-trail minimum, class / behind j 

V 

Approach speed 

D 

Common path length 


Arrival ROT on runway r 

&RAn 

Standard deviation of arrival ROT on runway r 

Rch 

Departure ROT 

GRDi 

Standard deviation of departure ROT 

Dor 

Distance to turn on departure on runway r 

v a 

Departure speed 

&VDi 

Standard deviation of departure speed 

Oxi 

Standard deviation of position uncertainty 

Ovi 

Standard deviation of approach speed 

G w 

Standard deviation of wind speed 

C 

Mean communications delay 

G c 

Standard deviation of communications delay 


Note: Subscripts indicate variation with aircraft class. ROT means runway occu- 
pancy time. 


Another rule requires that, if an arriving flight cannot make visual contact with the 
runway within 2 miles of the threshold, no departure can be released from that 
runway from the time the arriving flight is 2 miles out until it lands. This restric- 
tion also applies to a departure on a parallel runway within 2,500 feet. Advanced 
technologies may enable a departure release with confidence in this situation; if 
so, the user may either toggle this rule off or adjust the release distance to a value 
less than 2 miles. 

From the controller-applied separation, the module determines the time between 
arrivals of different aircraft classes. It computes a weighted average of these inter- 
arrival times to estimate the runway capacity. The module considers different op- 
erating strategies for a single runway (maximize arrival rate, maximize departure 
rate, balance arrival and departure rates) to generate a Pareto frontier of single 
runway capacity. 

Appendix A, the ASAC Airport Capacity Model and Delay Models Users Guide, 
includes sample input and output views. Appendix B details the derivation of the 
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Runway Capacity Module and how the inputs are translated into controller- 
applied separations. 

Multiple Runway Modules 

In some cases, the airport geometry is such that different runways can be operated 
independently in certain weather; in these cases, the capacity of the configuration 
can be determined from the capacities of the individual runways in that configu- 
ration. The technique for combining a set of runway Pareto frontiers to arrive at 
the configuration capacity curve is described in Appendix C. 

In other cases, operational interactions among different runways reduce the ca- 
pacity of the runways. In some cases this is modeled at the airport model level, 
using site- specific rules for determining the capacity of a configuration from the 
individual runway capacities. These site-specific rules are based on consultation 
with the air traffic controllers who work at the airport. Details are included in the 
descriptions of individual airport models; for example, see Appendix D for the 
description of the DTW capacity model. 

Other cases of runway interaction are common enough to warrant additional ana- 
lytical modules to compute the runway group capacity directly from the parame- 
ters that determine controller and aircraft behavior, as in the Runway Capacity 
Module. Two such modules are described in the following subsections. 

Departures on Closely Spaced, Parallel Runways 

The operational rules for parallel runways depend upon the centerline separation 
of the runways. If the runways are separated by no more than 2,500 feet (but at 
least 700 feet), a particular set of operational rules are stipulated by the FAA. 
Following convention, parallel runways separated by no more than 2,500 feet are 
referred to in this document as “closely spaced, parallel runways.” 

One restriction on closely spaced, parallel runways applies if both runways are 
being used to conduct departure operations. After a heavy jet departs, no depar- 
tures are permitted from either runway for 2 minutes. The user can instruct the 
ASAC Airport Capacity Model to take this rule into account, disregard it, or use 
the rule with a delay value of less than 2 minutes, as appropriate to the technology 
under consideration. 

Further restrictions on closely spaced, parallel runways are in place when one 
runway is used for departures, the other used for arrivals, and weather conditions 
do not allow visual contact by arrivals when they are 2 miles from the threshold. 

In this case departures must be held whenever an arriving aircraft is within 2 miles 
of the threshold. A special module is not needed for this case because the pair can 
be modeled as a single runway, with departures experiencing no runway occu- 
pancy time for arrivals. 
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The ASAC Airport Capacity Model 


A final set of restrictions on closely spaced, parallel runways occur when both 
runways are used for departures and departing aircraft are not visible 1 mile out. 

In this case separation requirements exist both for aircraft departing the same 
runway and for aircraft departing the other runway. This is referred to as 
“staggered departures.” The modeling technique used to derive capacity for this 
situation is described in Appendix E. 

Staggered Arrivals 

Under certain conditions, FAA procedures permit independent arrivals on parallel 
runways with centerlines separated by at least some minimum distance (currently 
700 feet). Such operations are referred to as “simultaneous.” As ceiling or visibil- 
ity degrades, safety considerations may require that more restrictive operational 
rules be enforced, which results in spacing requirements between the arrivals on 
the parallel runways. The modeling technique used to derive capacity in this case 
is essentially the same as the one used for staggered departures, and it is described 
in Appendix E. 

Parameters and Data Sources 

In the ASAC Airport Capacity Model, the aircraft are designated as being in a 
particular aircraft class. Aircraft in the same class have the same average values 
for approach speed, departure speed, arrival runway occupancy time, and depar- 
ture runway occupancy time. The standard deviations of these quantities account 
for both randomness in performance and variations from the mean within the 
class. 

Each aircraft class is also identified as being heavy or not. The wake vortex hold- 
time rules (see description in The Runway Capacity Module section) apply to air- 
craft departing after, or flying through the flight path of, a member of the heavy 
classes. 

For some airports runway usage is constrained to aircraft with a noise profile be- 
low some threshold. For these airport models, each class must be identified as 
noise track restricted or not. 

The other data that depends on aircraft class are the separation matrixes. Each 
separation matrix gives the required terminal area separation (in nautical miles) 
for each class combination of leader and follower. The user must specify three 
such matrixes: one for visual approaches (VMC1), one for radar approaches with 
visual completion (VMC2), and one for instrument approaches (IMC). 

The user completes the traffic specification by entering, for each aircraft class, the 
fraction of total traffic that the class contributes. 
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The generality of the aircraft descriptions permits a wide variety of analyses to be 
performed. The user may define aircraft classes to suit his or her needs. 

Airport-specific data, such as common path length and distance-to-tum on de- 
parture, are provided by ASAC as defaults. The default values are based on dis- 
cussions with air traffic controllers at each airport. 

Communications delay is also provided as a default. Communication delay mod- 
els the time from when a departure clearance is issued to when takeoff roll begins. 

Another value provided as a default is the standard deviation of wind speed. This 
standard deviation is based on typical variations within an hour, or between dif- 
ferent strata. The long-run variation in wind speed is not appropriate here, as the 
ASAC Airport Delay Model receives the actual surface-wind data for each hour 
from a weather database. 

A final user input to the ASAC Airport Capacity Model applies only for those 
cases where staggered departures or staggered arrivals may be used. (See Multiple 
Runway Module for a description of these modules.) The capacity for these run- 
ways is estimated by bounding the capacity from below and above. The user may 
set an integer computational parameter that controls the bounding process. The 
larger the value of the parameter, the more effort is put into making the bounds 
close, at the expense of additional computational time. ASAC provides a default 
value that gives a good tradeoff between accuracy and time. The ASAC output 
enables users to see the value of the upper and lower bounds, so as to make de- 
terminations about increasing or decreasing the parameter value. A full descrip- 
tion of the technique is contained in Appendix E. 

The ASAC user may modify the default values described in this section. 

Limitations of Runway Modules 

The ASAC Runway Capacity Module assumes that traffic appears in a random 
fashion. In fact, the controllers have some flexibility in sequencing arrivals to 
achieve desirable ordering (i.e., a heavy jet behind a small jet, rather than the re- 
verse). The controllers do not have unlimited flexibility in this regard, however, 
so modeling this aspect of air traffic control is best addressed at an airport-specific 
level. Not accounting for arrival or departure sequencing will result in slightly 
pessimistic estimates of capacity. 

Limitations of Airport Capacity Models 

Airports and weather are highly complex systems; a model with high fidelity 
would require both extensive development effort and long run times. Such models 
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The ASAC Airport Capacity Model 


already exist; however, to provide rapid indication of the likely effects of ad- 
vanced technologies, the ASAC model developers consciously constructed the 
airport models to reflect the predominant features of the arrival and departure en- 
vironment, while ignoring those rare situations (runway closures due to mainte- 
nance, for example) that will not have a significant long-term impact. 

It is not uncommon for airport capacity to be limited by ground traffic congestion 
near the gates. The ASAC models do not include a surface-movement component, 
other than reductions in capacity needed for arrivals to taxi across active runways 
in some configurations. The models needed to analyze delays due to gate conges- 
tion may find the unimpeded capacity results of the ASAC Airport Capacity 
model useful as input. 
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Chapter 2 


The ASAC Airport Delay Model 


The ASAC Airport Delay Model focuses on arrival delays. These delays account 
for the bulk of the users’ costs. The delay model uses the results of the ASAC 
Airport Capacity Model as an input. If users have an alternative means of com- 
puting capacity, they can present a file in the proper format to the ASAC Airport 
Delay Model without first running the ASAC Airport Capacity Model. The re- 
quired format is airport specific and is described in the documentation for each 
airport model. See Appendix A for an example. 

Computing arrival delays requires two relatively independent components. The 
first component is airport specific: determining the specific arrival or departure 
capacity curve point of each configuration, given the current weather data 
(precipitation, wind direction, wind speed, ceiling, and visibility) and demand, 
and selecting the best configuration. The second component is the same for all 
airports: using the selected configuration capacity and demand data to arrive at 
delay estimates. 

In the ASAC Airport Delay Model, an aircraft is considered delayed if it is not 
able to land as soon as possible. The ASAC Airport Delay Model Queuing Mod- 
ule computes a statistical estimate of delay due to congestion given the arrival 
demand and airport capacity. Even in those cases where demand is less than ca- 
pacity, some delay will be reported. This is because the statistical estimates ac- 
count for the random nature of queuing systems, in which arrivals in close 
temporal proximity will generate a delay, even though the long-run average ca- 
pacity of the system exceeds demand. Of course, the more that capacity exceeds 
demand, the smaller the estimate of delay. In fact, because of the hub-and-spoke 
scheduling practices common among the major airlines, arrivals tend to be 
bunched even more tightly than randomness alone would suggest. As such, the 
delay estimates of the ASAC system are optimistic. This optimism should not im- 
pair the ASAC Airport Delay model’s ability to compare capacity enhancements 
because the nature of the error is consistent in all scenarios. 

In the absence of congestion, some flights will arrive sooner than scheduled. The 
ASAC Airport Delay Model does not reduce the delay metric to account for early 
arrivals; it only tallies the minutes later than expected. 
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Model Results Versus Reported 
Performance Metrics 

Both the FAA and the airlines keep records of delays due to congestion in the air 
traffic system. Neither data source is likely to agree with the arrival delays re- 
ported by the ASAC Airport Delay Model because of the methods used to record 
delays. 

The ASAC Airport Delay Model accounts for delay in arrival time due to conges- 
tion. This delay may be incurred by speed reductions or vectors assigned by air 
traffic controllers, as well as by the less common occurrence of airborne holding. 

The FAA’s delay reporting system, OPSNET, only assigns delay if airborne 
holding is used by controllers, and then only if the time spent holding exceeds a 
minimum threshold (currently 15 minutes). Additional flight time caused by speed 
reductions or vectoring are neither recorded by controllers nor reported in 
OPSNET. 

The airlines’ delay information, as recorded in databases such as the Airline 
Service Quality Performance (ASQP), is based on scheduled segment length as 
the nominal flight duration. To maintain on-time performance, and to make sure 
that connecting passengers have enough time to transfer at the destination, the 
scheduled segment lengths are computed by the airlines to account for anticipated 
delays above the unimpeded flight time. Thus the ASQP delay information only 
reflects delays beyond that anticipated, as contrasted with the ASAC delay metric, 
which does not anticipate any delay. 

Although the delay metrics computed by the ASAC system are not directly com- 
parable to existing reported delay data, for purposes of evaluating the effects on 
system capacity, the ASAC metric is preferable. It is not biased for existing con- 
gestion, nor does it ignore delays because of their small duration or the means by 
which they were incurred. 

Configuration Selection 

The ASAC Airport Delay Model takes as inputs hourly weather data. These data 
consist of precipitation, wind speed, wind direction, ceiling, and visibility. The 
model uses airport-specific rules to determine which runways can be used, given 
the weather of that hour. For a configuration to be usable, all of its runways must 
be usable. The weather data also determine the capacity curve that describes each 
configuration. The model uses the hourly demand to determine the operating point 
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The ASAC Airport Delay Model 


on each configuration’s capacity curve that has the same ratio of arrivals to de- 
partures as the demand. The configuration selected as the one that would be oper- 
ated is the one with the greatest total capacity. Output to the queuing module is 
the demand and capacity for that hour. 

Arrival Delay Queuing Module 

The sporadic nature of air traffic demand, combined with the fluctuations in air- 
port capacity due to weather, results in a system that only reaches a scenario of 
constant arrival rate and constant acceptance rate (steady state in queuing terms) 
during those uninteresting early-morning hours when no one is flying. During the 
periods of interest, the conventional closed form queuing formulas, which rely on 
steady state assumptions, are of little value in estimating delays. ASAC requires a 
more dynamic model. 

The ASAC queuing module is an implementation of a published algorithm. 1 It 
uses a simplifying approximation to generate the time-dependent mean number in 
the system for a nonstationary queuing system, such as an airport. Empirical re- 
sults indicate that the method is sound. 

Parameters and Data Sources 

The hourly demand data for each airport are based on the 1993 schedule of aircraft 
operations published in the Official Airline Guide (OAG). For most large U.S. 
airports, which are the airports modeled in ASAC, the amount of traffic due to 
nonscheduled aircraft (military flights and general aviation) is an insignificant 
proportion of the total traffic mix. In those cases where discussions with air traffic 
controllers indicate that nonscheduled traffic may be significant, the OAG data 
has been modified to account for that traffic. 

The ASAC Airport Delay Model uses the FAA’s Terminal Area Forecast (TAF) 
data to estimate the growth in air traffic demand from 1993 to each following year 
(up to 2015). The user may choose a TAF traffic-inflation factor to increase the 
demand by some multiple of the 1993 baseline. The user may also input a traffic- 
inflation factor that is independent of the TAF. 

The weather data the model provides to the user are from actual observations at 
the airport in question. An entire year of historical data is provided for each air- 
port. The year is selected to be typical for that airport, in that the number of hours 
in each meteorological condition is roughly the same as the long-run average 
number of hours per year in each meteorological condition. The user may choose 
any beginning and ending day within the range of weather data available. 


1 Michael H. Rothkopf and Shmuel S. Oren, “A Closure Approximation for the Nonstationary 
M/M/s Queue,” Management Science, 25:6, June 1979, pp. 522-534. 
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Limitations of Airport Delay Models 


The ASAC Airport Delay Model does not account for delays attributable to airline 
scheduling practices that compact demand into small portions of each hour. As 
this assumption is the same in all scenarios, it should not affect comparative 
analysis. 

The model also does not account for the flight cancellations and schedule disrup- 
tions that occur in severe weather. The ASAC Airport Delay Model assumes that 
airlines intend all their arriving aircraft to arrive during their scheduled hours. As 
airport capacities are generally greatly reduced in such weather conditions, the 
ASAC Airport Delay Model will indicate a great deal of delay. In practice, this 
delay may have occurred on the ground and not in the air; however, as the occur- 
rence of such weather is relatively rare, the lack of fidelity in this regard should 
not appreciably affect the usefulness of the results. 
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Chapter 3 

Conclusions 


The ASAC Airport Capacity Model uses a set of parameters describing the rules 
of the air traffic system and the performance characteristics of the aircraft using it 
to compute runway capacities. The model combines runway capacities to generate 
capacity curves, which describe the tradeoff space between arrivals and depar- 
tures, for each operational configuration at the airport and for each category of 
weather at the airport that affects operational rules. 

The generality of the parameters allows the user to model the capacity impacts of 
changes to the basic characteristics of the system. The air traffic controllers 
working at each modeled airport evaluated the baseline cases of the model, using 
current operational parameters, and concurred on the model’s accuracy for their 
airport. 

The ASAC Airport Delay Model uses the results of the ASAC Airport Capacity 
Model, together with actual weather observations and demand patterns, to deter- 
mine the airport configuration and its capacity for each hour of the sample days 
selected. The resulting demand and capacity values are delivered to a queuing 
module, which computes the expected delay at the airport under those conditions. 

The ASAC modelers designed the capacity and delay models to give rapid results 
while maintaining sufficient complexity to account for the major features affect- 
ing the air traffic system. The resulting system is a powerful tool for analysts to 
use for quick evaluation of the delay reductions at existing airports that would re- 
sult from changes in the system parameters. In conjunction with the ASAC eco- 
nomic models, the ASAC airport models will provide a user-oriented system that 
facilitates benefit calculations of advanced air traffic tools and technologies. 
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Appendix A 

AS AC Airport Capacity and Delay Model 
Users Guide 


Airport Capacity Model 

The ASAC Airport Capacity Model is located in the “Model Wizard” and may be 
accessed at the following address: 

http://www.asac.lmi.org/access/model-wizard/ 

It is also accessible through the ASAC home page by making the following 
choices, in sequence: 

♦ select “Quick Response System (QRS)” 

♦ select “Enter” 

♦ select “Model Server” 

♦ select “Model Wizard” 

However the Model Wizard is accessed, the ASAC Airport Capacity Model is 
accessed by: 

♦ selecting “ASAC Airport Capacity Model” from the main dialog box 

♦ clicking on “CONTINUE” 

To perform an integrated capacity-delay analysis, using the ASAC Executive As- 
sistant is the preferred method of access. From the ASAC home page, select: 

♦ “Executive Assistant” 

♦ “Enter” 

♦ “Air Traffic Management Analysis” 

♦ “CONTINUE” 

When access to the capacity model is through the Executive Assistant, the session 
name and capacity model output will automatically be transferred to the ASAC 
Airport Delay Model. 
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Upon entering the capacity model, a brief description of the model appears, in- 
cluding a list of its four user input components: (1) aircraft class definitions; (2) 
aircraft separation matrices; (3) operating environment; and (4) runway. 

On this screen, the user is prompted for a session name. The session name will be 
used to identify all data and output files created during this session. A default 
name is provided. After entering a session name (or accepting the default), click- 
ing on “CONTINUE” will take the user to the next screen and through the user 
input components in the order they are listed above, to allow the user to customize 
the capacity model. 

Aircraft Class Definitions Component 

The Airport Capacity Model Aircraft Class Definitions File Locator (Screen 1) is 
the first screen to be displayed after selecting a session name. It offers the user the 
following three choices for loading information into the capacity model: 



Screen 1 


♦ “Find Aircraft Class Definitions File on Server” enables the user to select a 
file that is available on the ASAC system. This file may be used as is, or it 
may be edited and saved. After selecting “Find Aircraft Class Definitions File 
on Server” and clicking on “CONTINUE”, a list of available files is displayed, 
Select a file from the list box, then click on “CONTINUE”. The user may 
next click on “VIEW/EDIT FILE” or “CONTINUE”. If “VIEW/EDIT FILE” 
is selected, a screen with editable data fields, similar to Screen 2, will be dis- 
played. Move between editable fields by using the [Tab] key on the keyboard 
or by clicking on the desired field. At this point the user may make changes to 
the information or continue to view all the different screens of information by 
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selecting ’’VIEW/EDIT NEXT AIRCRAFT CLASS” until all the related 
screens have been viewed and/or changed as needed. 
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If at any point the user decides that changes are unnecessary, even though 
some have been made, clicking on the “RESET” button will reproduce the 
original information on that screen. 





The browser’s Back button can be used to return to editing a previous screen. 
’’VIEW/EDIT NEXT AIRCRAFT CLASS” must be clicked to save new 
changes. 


At the end of the editing process, select “FINISHED VIEWING/EDITING 
FILE”. The next screen allows the user to assign a name to the file. A default 
name based on the session name is provided. After entering a file name, if de- 
sired, click on “SAVE CHANGES” for any changes made. The file will be 
saved with the given name in the user’s ASAC directory. Alternatively, the 
user may click on “CONTINUE” if there are no changes. 

♦ “Build New Aircraft Class Definitions File” enables the user to develop a file 
with customized information, (Screen 2). Once all the different classes with 
pertinent information have been created, click on “FINISHED BUILDING 
FILE” . The next screen allows the user to assign a name to the file. A default 
name based on the session name is provided. After entering a file name, if de- 
sired, click on “SAVE FILE”. The file will be saved with the given name in 
the user’s ASAC directory. 

There is a limit of seven classes that may be defined when building a model on- 
line. To create a model with more than seven classes, the “Upload Aircraft Defi- 
nitions File to Server” option must be selected. 

♦ “Upload Aircraft Definitions File to Server” enables the user to use a file from 
their own computer. Up to 27 classes may be defined in an uploaded file. 
Follow the instructions that are listed on the screen for this process. 

Uploaded input file format. The uploaded aircraft class definitions file must be- 
gin with an integer giving the number of classes to be defined. Each class defini- 
tion is then given on a separate line, with fields separated by the tab character. The 
user may add comment lines to the uploaded file by inserting the character # in the 
first column of the line; any such lines are ignored by ASAC. 

The fields, in order, are: 

♦ class name 

♦ heavy jet flag (1 if heavy, 0 if not) 

♦ approach speed 

♦ standard deviation of approach speed 

♦ IMC arrival runway occupancy time 

♦ standard deviation of IMC arrival runway occupancy time 
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♦ VMC arrival runway occupancy time 

♦ standard deviation of VMC arrival runway occupancy time 

♦ departure runway occupancy time 

♦ standard deviation of departure runway occupancy time 

♦ departure speed 

♦ standard deviation of departure speed 

♦ percentage of this class in airport traffic 

♦ positional uncertainty to ATC. 

The class name may be any name the user selects. 

Heavy jet flag indicates whether wake vortex hold rules apply for departures after 
this aircraft class. 

Approach speed is the average speed over the length of the common path. De- 
parture speed is the average speed over the distance until which departures may 
turn. For on-line entry these values are in knots. For uploaded files the values are 
in miles per hour. The acceptable range of values for approach speed is 100 to 300 
knots. The acceptable range for departure speed is 0 to 300 knots. Any values out- 
side these ranges will generate an error. 

Arrival runway occupancy time is the time from when an arrival crosses the run- 
way threshold until it turns off the runway. There is one such value defined for 
VMC and one for IMC. Departure runway occupancy time is the time from when 
departure roll starts until the aircraft is no longer over the runway. Both times are 
in minutes. Acceptable values for arrival runway occupancy time are 0 to 3 min- 
utes. Acceptable values for departure runway occupancy time is 0 to 4 minutes. 
Any values outside these ranges will generate an error. 

Positional uncertainty to ATC refers to the controllers’ uncertainty in actual air- 
craft position given the position reported by radar to the control tower. A value of 
0.25 nautical miles reflects current technology. In uploaded files, this value should 
be in statute miles. A negative value will generate an error message. 

Percentage of this class in airport traffic gives the fraction of the airports opera- 
tions that involve aircraft of this class. ASAC checks that the sum over all classes 
is 1.00. 
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Standard deviation, for each item for which it is available, refers to the combined 
uncertainty in an individual aircraft’s performance and the variations among dif- 
ferent aircraft within a class. In all fields where a standard deviation is required, 
the standard deviation must be less than 1/5 of the corresponding mean value for 
that item. 

Once the steps for one of the above choices for inputting aircraft class definitions 
has been completed, the model will automatically lead the user to the next input 
component. 

Aircraft Separation Matrices 


As in the previous component, the File Locator (Screen 3) is the first screen to 
appear, offering the user the following three similar choices: 



Screen 3 


• “Find Aircraft Separation Matrices File on Server” enables the user to select a 
file that is available on the ASAC system. This file may be used as is, or it 
may be edited and saved. After selecting “Find Aircraft Separation Matrices 
File on Server” and clicking on “CONTINUE”, a list of available files is dis- 
played. Select a file from the list box, then click on “CONTINUE”. The user 
may next click on “VIEW/EDIT FILE” or “CONTINUE”. If “VIEW/EDIT 
FILE” is selected, a screen with editable data fields, similar to Screen 4, will 
be displayed. At this point the user may make changes to the information or 
accept the information by clicking on “CONTINUE”. The next screen allows 
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the user to assign a name to the file. A default name based on the session 
name is provided. After entering a file name, if desired, click on “SAVE 
CHANGES” for any changes made. The file will be saved with the given 
name in the user’s ASAC directory. Alternatively, the user may click on 
“CONTINUE” if there are no changes. 


Builder 








♦ “Build New Aircraft Separation Matrices File” enables the user to develop 
a file with customized information, involving separation of aircraft, by 
class, during different meteorological conditions (Screen 4). After entering 
data in the fields, click on “CONTINUE”. The next screen allows the user 
to assign a name to the file. A default name based on the session name is 
provided. After entering a file name, if desired, click on “SAVE FILE” in 
order to save the new file. 

♦ “Upload Aircraft Separation Matrices File to Server” enables the user to 
use a file from their own computer. This option must be selected if there 
are more than seven classes. Follow the instructions that are listed on the 
screen for this process. 

Uploaded input file format. The uploaded file should have one separation matrix 
for each of the three meteorological conditions. Each separation matrix consists of 
one row (line) per aircraft class. The order of the rows corresponds to the order in 
which classes were defined. The row for each class has a value for every class, 
separated by tabs. The value is the terminal area separation required when the row 
class is following each of the others. For on-line entry this value is in nautical 
miles. For uploaded files the values are in statute miles. The acceptable range of 
values is 1.5 to 7.0 nautical miles. Any values outside of this range will generate 
an error. 

The user may add comment lines to the uploaded file. Each comment line begins 
with the # character in the first column. 

Once the steps for one of the above choices for inputting aircraft separation matri- 
ces has been completed, the model will automatically lead the user to the next in- 
put component. 
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Operating Environment 


As with the other input components, a File Locator screen (Screen 5) is the first to 
appear for this input component. The following four choices are available: 



Screen 5 


♦ “Use Default Operating Environment File” brings the user to a screen that 
offers the choices “VIEW/EDIT FILE” or “CONTINUE”. At this point the 
user may elect to either view the available information by clicking on 
“VIEW/EDIT FILE” (which displays a screen similar to Screen 6),and 
then edit it if changes are deemed necessary. The next screen allows the 
user to assign a name to the file. A default name based on the session 
name is provided. After entering a file name, if desired, click on “SAVE 
CHANGES” for any changes made. The file will be saved with the given 
name in the user’s ASAC directory. Alternatively, the user may click on 
“CONTINUE” if there are no changes. 
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Screen 6 


♦ “Find Operating Environment File on Server” enables the user to select a 
file that is available on the ASAC system. This file may be used as is, or it 
may be edited and saved. After selecting “Find Operating Environment 
File on Server” and clicking on “CONTINUE”, a list of available files is 
displayed. Select a file from the list box, then click on “CONTINUE”. 

The user may next click on “VIEW/EDIT FIT F” or “CONTINUE”. If 
“VIEW/EDIT FILE” is selected, a screen with editable data fields, similar 
to Screen 6, will be displayed. At this point the user may make changes to 
the information. The next screen allows the user to assign a name to the 
file. A default name based on the session name is provided. After entering 
a file name, if desired, click on “SAVE CHANGES” for any changes 
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made. The file will be saved with the given name in the user’s ASAC di- 
rectory. Alternatively, the user may click on “CONTINUE” if there are no 
changes. 

♦ “Build New Operating Environment File” enables the user to develop a 
file with customized information. A screen similar to Screen 6 is dis- 
played. Once changes have been made click on “CONTINUE”. The next 
screen allows the user to assign a name to the file. A default name based 
on the session name is provided. After entering a file name, if desired, 
click on “SAVE CHANGES” for any changes made. The file will be saved 
with the given name in the user’s ASAC directory. Alternatively, the user 
may click on “CONTINUE” if there are no changes. 

♦ “Upload Operating Environment File to Server” enables the user to use a 
file from their own computer. Follow the instructions that are listed on the 
screen for this process. 

Uploaded input file format. Each field is on a separate line. The fields, in order, 
are: 

♦ communications delay 

♦ communications delay standard deviation 

♦ common approach path length 

♦ minimum distance before turn 

♦ wind uncertainty (standard deviation of wind speed) 

♦ wake vortex hold time 

♦ wake vortex after heavy departure rule flag 

♦ IMC Departure release distance 

♦ use IMC Departure release distance rule 

♦ meteorological condition definitions 

Communications delay is the time from when a departure clearance is issued by 
the air traffic controller until takeoff roll begins. Time is in minutes. Acceptable 
values are 0 to 1 minutes. The standard deviation may not exceed 1/5 of the mean. 
Any values outside these ranges generate an error. 
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Common approach path length is the distance over which aircraft fly in-trail in 
approaching a runway. Minimum distance before turn is the distance over which 
departures are in-trail from a runway. Units are in nautical miles for online entry 
and statute miles for uploads. Negative values generate an error. 

Wind uncertainty refers to the variability of wind speed over the common ap- 
proach path. It does not refer to variations in wind speed over time. The ASAC 
Airport Delay Model uses actual wind data; hence, wind variations over time (as 
well as temporal correlation) are accounted for directly. 

Wake vortex hold time is the time from when an aircraft that is “heavy” (see air- 
craft class definitions) departs until the time when the next departure clearance on 
the same runway may be issued. The time unit is minutes. The acceptable range is 
0 to 3 minutes. Any values outside these ranges generate an error. This value is 
used only if the wake vortex after heavy departure rule flag is set. For uploaded 
files, a value of 1 indicates that the hold time rule should be used; 0 indicates that 
this rule should be ignored. 

IMC departure release distance is the minimum distance from the runway 
threshold that an arriving flight may be while a departure is still permitted from 
that runway. This rule is only used when the runway is not visible to the arriving 
flight at two miles distance. Units are in nautical miles for online entry and statute 
miles for uploads. The acceptable range is 0 to 2 nautical miles. Any values out- 
side these ranges generate an error. This rule is used in capacity computations 
only if the use IMC departure release distance rule flag is set. For uploaded files 
a value of 1 indicates that the release distance rule should be used; 0 indicates that 
this rule should be ignored. 

The last data item only needs to be provided in uploaded files; the ASAC system 
adds this information to data files created online. The meteorological condition 
definitions give the ceiling and visibility minima for VMC1, VMC2, IMC1, and 
IMC2. The units are hundreds of feet for ceiling and hundredths of miles for visi- 
bility. The definitions are particular to each airport and should not be changed . For 
Detroit Metropolitan Wayne County Airport (DTW) the following line should be 
input, with values separated by the tab character: 

#Minimum ceiling visibility for VMC1, VMC2, IMC1, IMC2 
45 500 10 300 2 34 2 30 

The user may add comment lines to the uploaded file by inserting the character # 
in the first column of the line; any such lines are ignored by ASAC. 
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Runway Component 

This is the final input component of the ASAC Airport Capacity Model. The fol- 
lowing four choices are available at the first screen (Screen 7) within this compo- 
nent: 
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Screen 7 

♦ ’’Use Default Runway Configuration File” leads the user immediately to 
running the ASAC Airport Capacity Model (see Running the Airport Ca- 
pacity Model subsection). 

♦ “Find Runway Configuration File on Server” displays a screen that offers 
the user a choice of files. Select a file and click on “CONTINUE”. The 
next step will be to run the ASAC Airport Capacity Model (see Running 
the Airport Capacity Model subsection). 

♦ “Build New Runway Configuration File” enables the user to customize the 
model with information on minimum ILS ceiling and visibility (Screen 8). 
The user may assign a name to this new file by entering one in the box at 
the top of the screen. A default name based on the session name is pro- 
vided. Once changes have been made click on “CONTINUE”. The next 
step will be to run the ASAC Airport Capacity Model (see Running the 
Airport Capacity Model subsection). 
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♦ “Upload Runway Configuration File to Server” enables the user to incor- 
porate a file of information from their system into the model. Follow the 
instructions that are listed on the screen for this process. 
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Screen 8 - (Second half) 

Uploaded input file format. The default input file for DTW is reproduced below, 
along with the comments describing each field. All data values on the same line 
should be separated by the tab character. The underlined values in boldface must 
be entered exactly as shown in an uploaded file. Failure to do so will invalidate 
the delay model results. 

# number of runways 

8 

# the data for each runway appear on a single line as tab-separated val- 
ues 

# Runway name ILS minimum ceiling ILS minimum visibility Heading 

3L -1 11 34.3 
3C -1 11 34.3 
3R -1 11 34.3 
21L-1 11 214.3 
21C-1 11 214.3 
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21R-1 11 214.3 


27L2 50 274.7 

27R 2 45 274.3 

# Magnetic deviation 
- 5.6 

The ILS minimum ceiling and ILS minimum visibility for each runway give the 
values below which this runway cannot be used. Ceiling minima are in hundreds 
of feet; visibility minima are in hundredths of miles. A value of -1 indicates that 
this item does not constrain the runway usage. The units are hundreds of feet for 
ceiling, and hundredths of miles for visibility. 

Magnetic deviation gives the variation, in degrees, of local headings from true 
headings. A negative value is used for deviations to the West. 

Running the Airport Capacity Model 

Running the ASAC Airport Capacity Model is contingent upon completing all 
four input components of the model. The actual running of the model automati- 
cally follows the end of the fourth input component. 

The first screen to appear lists the files from the input components (aircraft class 
definitions, aircraft separation matrices, operating environment and runway) that 
will be used in the model. If any of these values are in error, the user can use the 
“Back” feature of their browser to return to a previous step, and correct the entry. 
Once the user accepts the files that will be used, click on “RUN AIRPORT 
CAPACITY MODEL”. 

The next screen (Screen 9) indicates that model execution has completed and lists 
the results of running the model, which are a file for the Airport Capacity Fron- 
tiers and an input file for the Airport Delay Model. Select the Airport Capacity 
Frontiers. Screen 10 gives a sample output. The number of configurations that are 
available will be indicated at the top, (in this case 4) followed by the runway con- 
figurations (example: 21L/21C/21R). 
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Each runway configuration has listed below it five sets of numbers, representing 
the capacity for that configuration in each of five weather conditions: YMC1, 
VMC2, IMC1, IMC2, and IMC3. The first number of each weather condition set 
is the number of points that can be found on the capacity frontier for that configu- 
ration in that weather condition. Each point on the capacity frontier is an (arrival 
rate, departure rate) pair. For example, the first number to appear immediately un- 
der the runway configuration 21L/21C/21R is a 5. This indicates there are 5 points 
for VMC1 on the frontier. The first point is “0,88”, indicating that when there are 
no arrivals, 88 departures per hour can be accommodated in this configuration. 

The second point is “39, 88”, indicating that 39 arrivals per hour can be accepted 
while maintaining the maximum departure rate of 88 per hour. In the first configu- 
ration (21L/21C/21R), the fourth and fifth set of numbers, for capacity in IMC2 
and IMC3, respectively, indicate 1 pair of points that are “0,0”. This indicates that 
this runway configuration is not usable in IMC2 or IMC3. 


Next, click on the “Back” button on your browser, to return to the screen that lists 
the output results of running the model. Select the Airport Capacity Model output 
file for the Delay Model. Screen 1 1 is an example of the results. 


location:] ftp: //ftp. asac. Imi. org/pub/Models/Outpu(/cnelson1 266. 1 309. din 



Screen 11 

This file will become the Input file for the delay model. For the DTW model, the 
delay model input is five sets of single runway capacity values, one for each me- 
teorological condition. Each set of values is tab delimited, and gives for that me- 
teorological condition the maximum departure rate; departure/arrival rate when 
arrivals are spaced to accommodate a departure between each pair of arrivals; the 
number of departures per hour that can be accommodated without increasing arri- 
val separations above the minimum required; and the maximum arrival rate. 

If the user wishes to provide their own capacity file to the ASAC Airport Delay 
Model, it must conform to the format described in the previous paragraph. 
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Airport Ph i . ay Model 

The user is automatically transferred to the ASAC Airport Delay Model if “Air 
Traffic Management Analysis” was selected from the ASAC Executive Assistant. 

The ASAC Airport Delay Model may also be accessed from the “Model Wizard” 
at the following address: 

http://www.asac.lmi.org/access/model-wizard/ 

or, through the ASAC home page by making these choices, in sequence: 

♦ select “Quick Response System (QRS)” 

♦ select “Enter” 

♦ select “Model Server” 

♦ select “Model Wizard” 

However the Model Wizard is accessed, the ASAC Airport Delay Model is 
accessed by: 

♦ selecting “ASAC Airport Delay Model” from the main dialog box 

♦ clicking on “CONTINUE” 

A very brief description of the delay model is given in the first screen. The user 
will also be asked for a “session name” which will then be assigned to the output 
data files created during this session. Clicking on “CONTINUE” at this point will 
start the user through the next few screens that develops the scenario under which 
the delay model is run. 

Develop Scenario 

The first several screens in the delay model request information or actions from 
the user in the order that they are listed below. As each process is complete, pro- 
ceed to the next screen by clicking on the “CONTINUE” button at the bottom of 
the screen. 

♦ the range of dates for the weather data; 

♦ the inflation rate to be used for air traffic projections (The default uses a 
schedule based on the flights listed in the 1993 Official Airlines Guide 
(OAG); traffic growth projections are from the FAA’s Terminal Area 
Forecast (TAF). The 1993 hourly demands from the OAG are increased 
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by the traffic growth percentage selected before input to the delay model. 
The actual demands used can be viewed as part of the delay model out- 
put.); 

♦ to either upload or find the AS AC Airport Capacity Model output file 
(using a tool called the Output File Locator. The format for uploaded files 
is described in the Running the Airport Capacity Model section); 

♦ to locate the Runway File for the airport (using a tool called the Runway 
File Locator). 

File Locator 

The File Locator procedure mentioned previously requires the user to make one 
from several choices for loading the runway file into the ASAC Airport Delay 
Model (Screen 12): 
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Screen 12 

♦ “Use Default Runway File” uses a file available on the system. The user 
will be able to see which files have been selected and to run the delay 
model in the next step by clicking on “RUN AIRPORT DELAY MODEL” 
(Screen 13). The model results will appear on the next screen. 
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Netsca pe - |Run th« ASAC Awpotl Delay Model] 


Run the ASAC Airport Belay Model 


Press Rim AIRPORT DELAY MODEL to run the model and create the output Se(s). 


Screen 13 


♦ “Find Runway File on Server” enables the user to access a file by using the 
Finder tool. This tool will list available files from which the user must se- 
lect one. Select one, then click on “CONTINUE” to proceed to running the 
delay model. Click on the button marked “RUN AIRPORT DELAY 
MODEL” (Screen 13). The model results will appear on the next screen. 

♦ “Build New Runway File” enables the user to customize the delay model 
inputs through use of a tool called Airport Delay Model Runway Configu- 
ration File Builder. The Builder will set up a file with ceiling, visibility 
and heading information by runway (similar to Screen 8). The user may 
assign a name to this new file by entering one in the box at the top of the 
screen. A default name based on the session name is provided. Clicking 
on “CONTINUE” will save these changes. The next step is to run the de- 
lay model by selecting “RUN AIRPORT DELAY MODEL” (Screen 13). 
The model results will appear on the next screen. 

♦ “Upload Runway File to Server” enables the user to incorporate a file from 
their own computer by selecting and following the instructions as they ap- 
pear on the screens. The format for uploaded files is described in the Air- 
port Capacity Model User Guide section. 
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Running the Airport Delay Model 

The ASAC Airpot Delay Model will generate two outputs; the user may select 
either by clicking the mouse on their choice: 

1 . Delay in Minutes 

2. Demand and Capacity Values used by queuing Engine 


Selecting “1” will display Screen 14: 
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Screen 15 

The data appears in the following format: 


Arrival Arrival Time Departure Departure Time 

Demand Capacity Interval Demand Capacity Interval 

Each line shows the values of demand and capacity used for one time interval, be- 
ginning with the first time interval of the day. The time interval is “1”, indicating 
that the time interval is 1 hour. 
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Appendix B 

ASAC Runway Capacity Module 


Overview 

We have developed detailed model of runway operations and capacities. The pa- 
rameters we used are listed in Table B-l. 


Table B-l. Key Airport Modeling Parameters 


Symbol 

Definition 

c 

Communications time delay 

Sc 

Variation in c 

D 

Length of common approach path 

D d 

Distance until departures may turn 

Pi 

Fraction of operating aircraft that are type / 

Rai 

Arrival runway occupancy time of ith aircraft 

8R a , 

Variation in R A , 

Rot 

Departure runway occupancy time of ith aircraft 

8Roj 

Variation in R a 

S 

Miles-in-trail separation minimum 

v, 

Approach speed of aircraft / 

8V, 

Variation in approach speed of aircraft / 

V Di 

Departure speed of aircraft / 

8 V D , 

Variation in departure speed of aircraft / 

SW, 

Wind variation experienced by aircraft / 

SX, 

Position uncertainty of aircraft / 

Pi a 

Time increment imposed by controller, / behind j 


We assume that each of the Sc, 8R Ai , SR Di , 8V t , SV Dj , <5W,, and 8X t variables are 
independent, normal, and random, with mean of zero and standard deviation of 
a RAO Ordo °vo a vDi °wo or °xi as appropriate. 
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Our analysis takes a “controller-based view” of operations. That is, we assume 
that a person controls the aircraft, introducing time (or, equivalently, space) in- 
crements in operations streams to meet all applicable rules (e.g., miles-in-trail re- 
quirements) with specified levels of confidence. For example, consider the 
arrival— arrival sequence of Figure B-l . 

Figure B-l. Time Phase for Arrivals When Follower 
Velocity > Leader Velocity 



Figure B-l shows the space-time trajectories of two arrivals. Zero distance is the 
beginning of the common approach path. In this model, the controller maneuvers 
the following aircraft so that it enters the common approach path a time p after the 
lead aircraft enters it. (The controller may actually achieve this by bringing the 
following aircraft onto the common path when the lead aircraft has advanced a 
specified distance along the path.) The controller chooses the time interval /j in 
light of his or her knowledge of typical approach speeds for the two aircraft and of 
disturbances affecting their relative positions — winds, position uncertainties, 
variations in pilot technique. The controller does this to ensure that miles-in-trail 
requirements and runway occupancy rules are met, with assigned levels of confi- 
dence. This action of the controller, together with information on statistics of air- 
craft operating parameters and the disturbances to arrival operations (such as 
winds and position uncertainties), leads directly to statistics of operations and of 
runway capacity. 

Arrivals Only 

We consider first the controller-based paradigm for arrivals only. Two cases are 
important. The first, illustrated by Figure B-l, occurs when the mean approach 
speed of the following aircraft exceeds that of the leader. 
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Follower Velocity > Leader Velocity 

For this case, the miles-in-trail constraint (distance) applies as the leader crosses 
the runway threshold. At that time, the leader’s position is D. We derive a condi- 
tion on the controller’s interval, p, to guarantee that the miles-in-trail requirement 
is met (i.e., that at the time the leader crosses the threshold, the follower is at least 
distance S away from the threshold, with a probability of 95 percent). 

The position of the lead aircraft at time t, X L (t), is given by 

X L (/) = dX L + (V L +8V l + SW L )t , [Eq. B-l] 

and the position of the following aircraft at time t, Xrft), by 

X F (0 = SX F + (V F + SV F + SW F ){t - p) . [Eq. B-2] 

The leader crosses the runway threshold at time tix>, given by 

D-5X l 

tw ~ V L +SV L +SW L ' 


At time t LO , the follower is at X F (t LO ), given by 


X F {t LO ) = SX F + (V f + SV F + ) 


D-SX, 


k v l +sv l +sw l ^ j 


[Eq. B-3] 


We wish to derive a condition on p, to make D-X F (t LO ) >S with a probabihty of 
at least 95 percent. To keep the problem tractable, we assume that all disturbances 
are of first order and make Equation B-3 linear. When made linear. Equation B-3 
becomes 


X F (t w ) = SX F + 


DV r 


V L 


1 + 


sv F +&v F sx L dv L +dw L 


V F 


D 


V L 




1 + 


Vf J 


. [Eq. B-4] 


In this linear approximation, X F (t LO ) is a normal random variable of mean 


DV r 


HV F 


and variance 


= 


D 2 V f 


r <? 2 VF 




WF 


V 


+ - 


XL 


D 


-h- 


VI. 


+ a 


WL 


+ p 2 V 2 


VF 


+ a 


WF 


+ c 2 xf . [Eq. B-5] 
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The condition that D-X F (t LO ) > S, with probability at least 95 percent, may then 
be stated as 


or 


— HV F + 1.65 w ^ D-S 
V L 


[Eq. B-6] 


D D-S 1.65(7, 

u > + L 

V V V 

V L V F V F 


[Eq. B-7] 


Equation B-7 gives, in essence, the desired condition. As that equation stands, p 
appears on both sides of the inequality. Straightforward manipulations lead to an 
explicit condition on fi, which may be written 


A + Va 2 B 2 + C 2 (l-fl 2 ) 
l-B 2 


[Eq. B-8] 


where 


A = 


_D 

V L 


D-S 


[Eq. B-9] 


B 2 



[Eq. B-10] 


and 


C 2 


_ 1-65 2 

D 2 V 2 

(a 1 _i_ /t 2 

° VF ° WF 

0 2 xl 

, & VL & WL 1 

+ct! ] 

r 

11 

Vl 

{ V? 

+ D- 

b 

+ ° XF 


[Eq. B-l 1] 


To determine numerical values of the smallest fl that meet Equation B-7, the it- 
erative scheme 


B„ +l 


D 1 -65(7, (/Q 


where <T, (jix) is defined by Equation B-5, is convenient. 
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Now, let us develop a condition on p that will guarantee that the follower aircraft 
does not cross the runway threshold until the leader has left the runway, with a 
probability of 98.7 percent. The leader will exit the runway at time 


1 lx ~ 1 lo + 

and the follower will cross the threshold at time tpo, given by 


^ FO 


V f +5V f +SW f 


+ p. 


[Eq. B-12] 


Linearizing as above, we find that in the linear approximation t FO -t LX is a normal 
random variable with mean 


_D 



-R 


AL 


and variance 


— 2 

° 2 ~ y 2 

V F 


D 2 


+ 


(J VF + O wf ^ 


+ - 


D 

vi 


( °\l 

D 2 




v! 


+ <y, 


[Eq. B-13] 


It follows that the condition on p for the follower not cross the threshold until the 
leader has exited the runway; that is, t FO -t LX > 0 with a probability of 98.7 per- 
cent, is 


D D 

p> + R 

V V 

v L v F 


al 2.215<t 2 . 


[Eq. B-14] 


The controller will, in effect, impose that value of time interval p that is the 
smallest p satisfying both Equations B-7 and B-14. 


Given p, the time between threshold crossings of successive arrivals is, in our ap- 
proximation, a normal random variable of mean 


_D 

V, 





[Eq.B-15] 


and variance 


a 


2 

3 


D 2 

f 2 2 2 

V XF . G VF^~°WF 

D 2 

_1 

( o\l °Il+° 2 wl) 

v 2 

V F 

D 2 V 2 

\ U V F ) 

vl 

{ D 2 V l J 


[Eq. B-16] 
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Follower Velocity < Leader Velocity 


When the follower’s approach speed is slower than the leader’s, in the controller- 
based view the controller will bring the follower onto the common path after the 
leader has advanced a distance S along it, as illustrated in Figure B-2. 

Figure B-2. Time Phase of Arrivals when Follower 
Velocity < Leader Velocity 
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In this case, the positions of the two aircraft as functions of time are again given 
by Equations B-l and B-2. The miles-in-trail requirement is now that X L (ji )- 
X F (ji) > S, with a probability of at least 95 percent. As 

X L (/i) - X F (ji) = 8X l + (V l +SV L +SW L )fl- SX F [Eq. B-17] 

is a normal random variable of mean V L /j and variance 

(O yl + & wl) "F & XF + G xl » [Eq. B-18] 

it follows that the condition that the miles-in-trail requirement is met with a 
95 percent confidence, is 

H >A + 1.65^-. [Eq.B-19] 

*/. V L 
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Equation B-19 may be written as a single condition on p using Equation B-8 by 
replacing Equations B-9, B-10, and B-l 1 with the following new definitions: 



B 2 = 1.65 2 


a VL +a 


WL 


and 


C 2 = i6 5 2 £k±^k. 

vl 

The condition that the single-occupant rule is met with 98.7 percent confidence is 
derived exactly as is that condition for V F > V L (i.e., the condition expressed in 
Equation B-14). Also, in the present case, Equations B-15 and B-16 yield the 
mean and standard deviation of interarrival times, given p. 

Arrival-Departure-Arrival-Departure Sequences 

We can readily translate the results of the previous subsection to results for opera- 
tions where there is a departure between successive arrivals by replacing Ra l with 
Ra l + Rdf, where Rdf is the departure runway occupancy time of the interven- 
ing departure aircraft. This case is illustrated by Figure B-3. 

Figure B-3. Time Phase of Arrivals with Intervening Departure 
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It may be desirable to consider the effect of a communications lag c on the depar- 
ture. If so, then Ra l is replaced by Ra l + c + Rdf- 

Statistics of Multiple Operations 

At this point, we have expressions for the means and variances of normal random 
variables representing interarrival times for a variety of cases. Now we will use 
these to generate statistics of multiple arrivals, or multiple arrivals and departures, 
that we can then use to produce capacity curves for single runways and combina- 
tions of runways. 

First, we consider the statistics of sequences of arrivals only. Statistics of the 
overall interarrival time is determined by the mix of aircraft using the runway, 
with their individual values for the aircraft parameters of Table B-l. The fraction 
of the aircraft of type i in the mix is p r The previous results give the interarrival 
time for each pair as a normal random variable. Let the mean and standard devia- 
tion for aircraft type i following aircraft type j be and cr i; -, respectively. Then 
the distribution function for overall interarrival time is 


PAA(t) = 'Z'Lp i PjN(f,H l] ,a ij ), [Eq. B-20] 

* j 

where N(t\ jl, a) denotes the normal probability distribution function. Obviously, 
the distribution of interarrival times is not normal. An example of an interarrival 
time distribution of the Equation B-20 type is shown in Figure B-4. 

Figure B-4. Interarrival Time ( Minutes ) 
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As Figure B— 4 shows, the interarrival time distribution is not necessarily mono- 
modal. 

The mean and variance of the interarrival time distribution. Equation B-20, can be 
computed straightforwardly. The results are 

< > = £ X PiPjfy [Eq- 62 1 1 

I i 

and 

var (t M ) = '2j£p,p J (<ri +^)~<t AA > 2 [Eq. B-22] 

i j 

In principle, the analyst can compute exactly the distributions of total arrival times 
for all arrival sequences of arbitrary length and find exact values for the number of 
arrivals that can, with assigned confidence, be accommodated in 1 hour. These 
calculations involve sums of many terms, however, and this motivates a search for 
useful approximations. 

Sums of normal random variables are normally distributed, and it is tempting to 
approximate the distribution of sequences of many arrivals in such a way. An all- 
arrival sequence with J cases of aircraft of type i following aircraft of type j has a 
normal distribution whose parameters are easily computed. If the J were chosen 
so that Jjj = PjPjM, where M is the sum of the J tj , the resulting normal distribu- 
tion would be a good approximation for the distribution of long arrival sequences. 
Unfortunately, for the aircraft mixes at some airports, some of the p i values are 
only a few hundredths, so M would have to be several thousand for this approxi- 
mation to be accurate. 

Nevertheless, because much of our work to this point has been approximate, con- 
sidering this “very large sequence”-limiting case does not seem unreasonable. In 
this approximation, then, the time t M of M interarrival times has a normal distri- 
bution of mean M<t AA > and variance 


var(/ M ) = J PiPjtf ■ [Eq- B ' 23 ] 

i j 

This result suggests approximating the distribution of interarrival time with a 
normal distribution of mean <t AA > and variance v , given by 

[Eq.B-24] 
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The approximation of Equation B-24 can be used to compute the number of arri- 
vals that can be accommodated in 1 hour with 95 percent confidence. That num- 
ber M is determined by the condition 

(M* - 1) < t M > +1.65^/(M* — l)v, < 60 

(only M *- 1 interarrival times are required for M* arrivals), 1 which leads to the 
all-arrival capacity of a single runway as M* = w 2 + 1, where w is given by 


w = - 


1.65^ If 1 . 65 ^ 


2 <^ > 


V 2< ^ >) 


+ - 


60 


<‘aa > 


To compute the expected number of arrivals we use 


< M > 


60 


<t A 


Runway Capacity Curve 

At this point, we have one point on the single runway capacity curve, the one cor- 
responding to all arrivals and no departures. We can generate others. 

The distribution function of Figure B-4 suggests that there is a significant prob- 
ability of interarrival times being large enough to accommodate a departure. We 
can compute the number of “free” departures in a stream of < M > arrivals. 
These are departures that can be accommodated without introducing additional 
separations into the arrival stream, in this way: the distribution of interarrival time 
is given by Equation B-20. We assume that departure runway occupancy time, 
arrival runway occupancy time, and communications delay are normal random 
variables of means <R D >, <R A >, <c> and standard deviations Gr D , o ra , and o c . 
The difference t-R D -R A -c , where t is the interarrival time, is the excess time 
when a departure is released between two arrivals. (The lead arrival exits the run- 
way in time R A , and the departure exits the runway in time Rp+ c.) The distribu- 
tion of the excess time is 


p(t-R A -R D -c)=N(-<R A >-<R D >-<c >,Jo~^+a^+af ) , or ) 

* j 

=SS m N(H-<R D >-<R A >-<c>,{c 4 + ^+oiL+q 2 ), [Eq.B-25] 

*' j 


lr This derivation is valid for a single hour, considered in isolation. To compute the long-run 
average replace M* - 1 by Af*. 
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AS AC Runway Capacity Module 


where the symbol <8> denotes convolution. 2 The probability that t-R D -R A -c is 
positive, i.e., the probability that a free departure is possible between any random 
pair of arrivals, is given by 

P + = 1 -llPtPjCiO,*- < R d > - < R a > - < c >, 

i j 

1 / °rd +<7 ra +<7 c +a l) [Eq. B-26] 

where C{t, (i, o) is the normal cumulative probability function. This value is 
readily computed. Then the analyst may determine the number < N > of positive 
values t-R D —R A -c to be expected, in < M > draws, from the binomial distribu- 
tion for probability p+. 

Under IMC2 or IMC3 weather conditions, current FAA procedures require that 
departures be held if an arriving aircraft is within a certain distance of the runway 
threshold. (This distance is now 2 miles.) In our model, this has the effect of re- 
ducing the time available for “free” departing aircraft. Since the trailing arrival 
travels a distance less than the full length of the common path, the uncertainties 
embodied in Equation B-26 are also reduced. 

The appropriate modification to Equation B-26 is to reduce the mean time avail- 
able for free departures by 


Dj_ 

y F 


where Dj is the distance from threshold after which departures must be held. The 
variance in Equation B-26 is reduced by 


D t (2x D- D t ) 


o; F + <r. 


2 "\ 


WF 


A third point on the capacity curve, the point of equal numbers of arrivals and de- 
partures, may be computed by considering sequences of repeated arrival-departure 

: To account for variations in departure runway occupancy time (ROT), the analyst may re- 
place the single normal distribution of ROT with the distribution of departure ROT that would be 
found with K classes of departing aircraft, each with its own normal distribution of departure ROT. 
That is, 

K 

^ Di) ' 

1 

where q, denotes the fraction of departing aircraft that are of type i. 
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pairs, as described in the section Arrival-Departure-Arrival-Departure Se- 
quences. 

While the expressions for interarrival times and runway capacities developed pre- 
viously in this appendix are somewhat lengthy, they are readily evaluated numeri- 
cally. 

Departures 

Considerations similar to those for all arrival sequences, as treated in the section 
Arrivals Only, must be taken into account to develop to statistics of departures. 
The basic departure situation is shown in Figure B-5. 

Figure B-5. Time Phase of Departures 


6 r 



We model the trajectory of a departing aircraft by specifying its position, x{t), in 
terms of the parameters V D and R D , in this way: 


*(0 = 1 


1 V 

--P-t\0<t<R 

2 R d 


V D t--V D R D ,t>R D 


This model approximates an actual takeoff roll and climb out by a trajectory with 
constant acceleration from rest to departure speed V D , occurring in time R a , fol- 
lowed by continuing departure at constant speed V D . 


B-12 



ASAC Runway Capacity Module 


We model controllers’ actions on departures by the interdeparture time interval p, 
which is the time interval between the start of the lead aircraft’s takeoff roll to 
when a departure clearance is issued to the following aircraft. (The following air- 
craft begins takeoff role at time p + c, where c models the delay to move into po- 
sition.) We assume that, in effect, controllers adjust p to give specified confidence 
that miles-in-trail requirements, and other separation requirements, are met. 


Here again, the required control input varies, depending on whether the following 
aircraft is faster or slower than the lead aircraft. In the case of a faster follower, 
the constraining condition is that the separation requirement be met as the lead 
aircraft exits the system. At that time, the displacement of the lead aircraft is D d , 
the distance to turn on departure. The displacement of the following aircraft must 
not be greater than D d -S d , where S D is the minimum separation. After lengthy 
but straightforward steps, the controller finds that meeting this condition with 
95 percent confidence imposes the condition 


D d 1 / \ 

A* - 77~ + r ( R DL - R df ) 


Dd ~ S D 1.65 i 

- — c + v var 


r DL 


V, 


DF 


’ DF 


on p. The quantity var in the inequality just above is 


var = 


Dp . \ ( n _ p \ 1 I 2 , 2 \ , ^PF n 2 ®VDL + 

\j 2 ' R DF/ H C \ ^VDF + y 2 Dp y 2 


v v DL 

+ ^A/prGfiDF + VpF®C • 


V 

I V _DF 2 
a ® RDL 


The inequality may be reduced to an equivalent, exphcit condition on p. For nu- 
merical work, we find that iterative methods give the required values of p con- 
veniently. 


When the follower departs more slowly than the leader, the separation minimums 
apply as the follower lifts off, unless D d is sufficiently short that the leader can 
exit the system before the follower completes the takeoff roll. Applying the sepa- 
ration minimum as the follower lifts off leads to the condition 


Rr 


p> 


- + 


1 V, 


DF 


2 V r , 


-1 


\R r 


■c + - 


■ + 


1.65V var 1 (inequality A) 
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where 


var 1 — /i -h c + R df R dl 
v l 


\ 2 <&L 


+ ( 7 . 


WL 


■ + 


r DF 


2V, 




DL / 


1 

+ 4 


7 rr 2 i - 2 \ 

2 r>2 ^VDF 

^KDL **DF r / 2 


V 


r DL 


+ CT, 


c • 


7 


Alternatively, the controller might impose a value of /j that caused the follower to 
lift off just as the leader exited the system. That would lead to 


D r 


D 1 y D 2 d a VDL 

M — y h ^ Rdl ~ R df ~ c ^ 1.65^ ~ 


+ <T. 


WL 


1 


+^c 


DL 


DL 


(inequality B). 

Controllers would impose the less restrictive of inequality A or inequality B. Fi- 
nally, the single-occupant rule must be respected, which leads to 

fl > R dl - c + 2.215 + a c (inequality C). 

For our model, when the follower is slower than the leader, we choose 

B - max [min^, Bb), Be], 

where /i i is the lower bound on fl resulting from inequality i. 
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Appendix C 

Combining Pareto Frontiers 


The ASAC Airport Capacity Model develops capacity Pareto frontiers for indi- 
vidual runways and groups of interacting runways. The Pareto frontier describing 
the overall capacity of a configuration is derived from the combination of the 
component runway and runway group frontiers active in that configuration. In this 
appendix, we describe a general methodology for combining Pareto frontiers for 
configuration capacity. 

Pareto frontiers can be combined pairwise, to arrive at a capacity frontier for the 
combination. In the event that a configuration has three runways or runway 
groups, the overall capacity Pareto frontier is arrived at by first combining any two 
Pareto frontiers, then combining the composite frontier just obtained with the 
third runway group’s frontier. For more than three runway groups, each additional 
runway group frontier is combined with the composite frontier of the previously 
treated groups, until all have been incorporated; at each step, a pairwise combina- 
tion is made. 

Combining a pair of capacity Pareto frontiers is the fundamental operation in ar- 
riving at a configuration capacity Pareto frontier and is the subject of the remain- 
der of this appendix. 

A Pareto frontier is described by a set of points in the arrival-rate/departure-rate 
plane. We assume that the frontier is piecewise linear. The combined frontiers will 
also be piecewise linear. Let (a,, d ( ); i = l...m and (bj,ej); j=l...n be the set of 
(arrival, demand) rates of the Pareto frontiers to be combined. For each point we 
can compute the right-hand slope of the frontiers s i and f ., respectively: 

d M ~ d i 


For the last pair on the frontier, the right-hand slope is undefined. A special value 
is stored in s m to indicate this. The values t- are computed analogously. The ab- 
solute value of the slope indicates the rate at which the controller must decrease 
departure capacity to increase arrival capacity from that operating point. 

We begin constructing the combined frontier by generating the linear segments 
that represent the operational tradeoffs allowed by the constituent capacity curves. 
The upper envelope of these segments is the combined Pareto frontier. Each po- 
tential operational segment begins with a point (c k ,f k ),f k = d i + e } and 
c k = a j + bj for some pair (ij). When operating at this point, the controllers have 
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two options for increasing arrival capacity at the expense of departure capacity: 
increase arrival rate from a, on one complex or increase arrival rate from bj on the 
other complex. The best choice is the former if s ( > t- and the latter otherwise. 

The end point of the segment is chosen to reflect this best possible tradeoff deci- 
sion, i.e., a i+] bj,d i+l + e } in the former case and a i + b j+i ,d i + e J + ] in the 
latter. We also store the slope of the segment in its data structure, to avoid the 
need for recomputing it. We will denote this slope u k . A pointer to the segment 
beginning at this segment’s endpoint is also stored. (If there is more than one such 
segment, the one with the largest slope is pointed to.) We compute potential op- 
erational segments for each possible pairing of frontier points from the constituent 
capacity curves; thus the number of such segments is N = mxn. We sort the seg- 
ments by the left-hand edge arrival rate for later convenience. 

Determination of the upper envelope of operational segments, the combined 
Pareto frontier, begins with the leftmost segment, the one beginning at 
/, = d x + e x and c, =a x + b x . This point typically represents the all-departure ca- 
pacity of the complexes being combined. For convenience we will denote this 
segment’s endpoint by (c 2 ,/ 2 ). This segment is the candidate segment for the 
combined Pareto frontier. We then examine all other segments 

J that have arrival rate ranges overlapping the candidate seg- 
ment. These are easily identified, as their left-hand edge arrival rate c k is less than 
the right-hand edge arrival rate c 2 of the candidate segment. One of four cases can 
occur: 

1 . The right-hand end of the examined segment has an arrival rate c k , that is 
less than right-hand end of candidate segment, and the departure rate of 
examined segment’s right-hand edge /*■, which is less than the departure 
rate of the candidate segment for that arrival rate d\ + u\{ck—c\). In this 
case, the examined segment is dominated: it lies entirely below the candi- 
date segment. The examined segment can be discarded from any further 
consideration. 

2. The right-hand end of the examined segment has an arrival rate c^, that is 
less than the right-hand end of the candidate segment, and the departure 
rate of the examined segment’s right-hand edge /*■, that is greater than the 
departure rate of the candidate segment for that arrival rate, 

d\ + Hi(c*— ci). In this case, the examined segment crosses the candidate 
segment. This is the most complex case, and we proceed as follows: 

♦ Compute the point of intersection between the two segments. 

♦ Set the left-hand edge of the examined segment to be this intersection 
point. (The portion of the examined segment to the left of the intersec- 
tion is dominated and can be ignored.) 
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Combining Pareto Frontiers 


♦ Set the right-hand edge of the candidate segment to be the intersection 
point, and set the candidate segment’s pointer to the (modified) exam- 
ined segment. 

♦ Create a new segment from the intersection point to the previous end- 
point of the candidate segment. The slope of this new segment is the 
same as that of the candidate segment, and its pointer is set to the old 
pointer value of the candidate segment. This new segment, as it ex- 
tends beyond the examined segment, may still be part of the combined 
Pareto frontier. 

We continue comparing the revised candidate segment to unexamined 
segments that overlap its new extent. 

3. The right-hand end of the examined segment has an arrival rate c* , that is 
greater than the right-hand end of candidate segment, and the departure 
rate of the examined segment at the right-hand edge of the candidate seg- 
ment c 2 , that is less than the departure rate of the candidate segment for 
that arrival rate dy + u k (c 2 -c k )< f 2 . In this case, the overlapping portion of 
the examined segment is dominated, but the remaining extent of the ex- 
amined segment may still be part of the Pareto frontier. We do nothing in 
this case. 

4. The right-hand end of the examined segment has an arrival rate c k , that is 
greater than the right-hand end of the candidate segment, and the departure 
rate of the examined segment at the right-hand edge of the candidate seg- 
ment c 2 , is greater than the departure rate of the candidate segment for that 
arrival rate dy + u k (f 2 -c k )< / 2 . In this case, as in case 2, the examined seg- 
ment intersects the candidate segment. The difference is that the portion of 
the candidate segment to the right of the intersection is dominated and can 
be discarded. The procedure is the same as in case 2, except that the last 
step is not performed. 

Once we have examined all overlapping segments, we know no segments are 
above the candidate segment. The endpoints of the candidate segment (perhaps 
modified during the search) are on the combined Pareto frontier. The next candi- 
date segment is the one pointed to by the old candidate segment. The procedure is 
repeated for this new candidate segment. The procedure will end with a segment 
whose endpoint is the arrival capacity of the combination. After examining this 
segment, all other segments will either be part of the Pareto frontier or will have 
been discarded as dominated. 

When the form of the Pareto frontier for the complexes to be combined is known, 
it is possible to specialize and simplify the procedure described here. Such spe- 
cializations are used where possible. 
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Appendix D 

Detroit’s Airport Capacity Model 


FAA controllers and managers at Detroit Metropolitan Wayne County Airport 
(DTW) informed us that the airport almost always operates in one of four configu- 
rations: 21L/21C/21R, 3L/3C/3R, 27L/27R, or 27L/27R/21R. (Operations in the 
last configuration require traffic on runway 27R to hold short of runway 21R.) 
Prevailing winds cause the 21 configuration to be by far the one most often used. 

In the 21 configuration, 21L is used exclusively for arrivals and 21C exclusively 
for departures. The arrival and departure mix on 21R is adjusted to best accom- 
modate demand. The 3 configuration is operated similarly, with 3R devoted to 
arrivals and 3C devoted to departures, with the mix on 3L adjusted in view of de- 
mand. 

When the 27 configuration is in use, 27L is devoted to arrivals and the mix on 
27R is adjusted to best accommodate demand. For the 27L/27R/21R configura- 
tion, both 27 runways accommodate arrivals only; runway 21 R is devoted to de- 
partures. This configuration is used only in dry weather. 

In most cases where a runway is used to balance capacity, the achieved arrival and 
departure rates on that runway are determined from the runway’s Pareto frontier in 
such a way that the proportion of arrival demand and departure demand presented 
is equal to the achieved proportion of arrival and departure demand. 

The exception to the above rule is when the 3 configuration is used. In the 3 con- 
figuration, two runways are used for arrivals (in IMC1 or better) and one for de- 
partures when the balancing mix would require an arrival rate exceeding the 
departure rate by more than 20 percent. When arrival demand is not so heavy, this 
runway is managed for an equal number of arrivals and departures. This scheme 
reflects the practice of controllers at DTW. 

For an example of how the model determines capacity as a function of weather 
and demand, consider the following situation. The ceiling is 600 feet and visibility 
is 5 miles. Due to the low ceiling, we are in IMC1. Winds are from the north 
(blowing from heading 340°) at 13 knots. This creates a tailwind on the 21s in ex- 
cess of the maximum of 8 knots, so the configuration 21L/21C/21R is not usable. 
During a departure push the arrival demand during the next hour is 23, and the 
departure demand is 64. The arrival/departure capacity of the usable configura- 
tions is 3L/3C/3R: 64/78; 27L/27R: 33/47; 27L/27R/21R: 66/47. Since the 
3L/3C/3R configuration is best able to meet the predominant demand, departure 
traffic, it is the one selected. 
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Now consider the same weather conditions when there is an arrival push. Arrival 
demand in the next hour is 67, and departure demand is 25. During an arrival 
push, two runways of 3L/3C/3R are dedicated to arrivals, rather than using one for 
arrivals and one for both arrivals and departures as in a departure push. In the 
27L/27R configuration, the runway used to balance traffic between arrivals and 
departures will be skewed toward accommodating arrivals. The configuration ca- 
pacities are 3L/3C/3R: 66/47; 27L/27R: 66/24; 27L/27R/21R: 66/47. We see that 
even for the same weather conditions, the demand affects the configuration ca- 
pacity. All of the configurations give equal arrival capacity during this arrival 
push. The 3L/3C/3R and 27L/27R/21R configurations allow the greatest departure 
capacity, so we would choose one of those. Since they have equal capacities, the 
choice is arbitrary. 



Appendix E 

Staggered Departure and Arrival Models 


In this appendix we describe the ASAC Airport Capacity Model algorithm used to 
estimate the capacity of a parallel runway pair when there are spacing require- 
ments between both aircraft using the same runway and between aircraft using one 
runway and aircraft using a parallel runway. This can occur when both runways 
are used for departures or when both runways are used for arrivals. 

Unlike separation requirements for single runways, separation requirements in this 
situation between aircraft approaching the same runway cannot be derived by ex- 
amining aircraft class pairs in isolation; the interdependence of traffic on the two 
runways requires, in general, knowledge of the entire sequence of operations to 
determine the separation required between any two aircraft approaching the same 
runway. 

Since exact separations cannot be determined, except for a specific sequence of 
operations, the algorithm constructs upper and lower bounds on the separation 
time required between successive operations on one runway of the pair. The 
bounds are computed for each combination of following aircraft class and leader 
aircraft class (as in the single runway model). The bounds take into account the 
interaction with traffic on the other runway. 

A user-controllable parameter determines how many historical operations are con- 
sidered, and thus how much refinement is put into determining the separation 
bounds, so that capacity can be estimated to any desired degree of precision (at the 
expense of additional computation time). The capacity bounds of the runway are 
computed on the basis of the weighted average time between operations; the 
weighting factors account for the traffic mix on the targeted runway. Since we as- 
sume that operations alternate between runways, the capacities of both the tar- 
geted runway and the “other” runway will be the same. We can exploit this 
symmetry by computing the capacity bounds twice, one using each runway as the 
target. The computed bounds will generally differ, leading us to identify a best 
lower bound and a best upper bound on estimated capacity. 

Here we discuss the capacity bounding algorithm from the perspective of depar- 
tures. The staggered-operations capacity algorithm for arrivals is completely 
analogous. 
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Modeling Departure Capacity of a Parallel 
Runway Pair 


In modeling the interdeparture times on the target runway, we assume that a de- 
parture has just occurred on the other runway. To capture the separation times re- 
quired between two aircraft on the target runway (aircraft of type i, following an 
aircraft of type j, which is next to depart), we need to consider also the aircraft of 
type /, which has just departed on the other runway, and the aircraft of type k, 
which is due to depart the other runway after the aircraft of type j departs the run- 
way under consideration. The departure sequence is l,j, k, i. For conciseness we 
will refer to an aircraft of type x as simply aircraft x. 

We define p(i,j,k,l) to be the average time separation (in minutes) that the con- 
troller will apply to aircraft i following aircraft j on the same runway, when air- 
craft l has just departed the other runway and aircraft k is next to depart the other 
runway. We compute both upper and lower bounds on this separation. 

The separation (in minutes) between i and j that we use to compute the runway’s 
capacity is the weighted average 

k,l 

where p xk (p xl ) is the probability of aircraft k(l) on the other runway. Upper 
(lower) bounds on p p (i,j) are computed using the upper (lower) bounds on 

p(i,j,k,l). 

The hourly runway capacities are estimated by 

60 

capacity = ^ , 


where />, and p, are the probability of i and j on the targeted runway. Lower 
(upper) bounds on capacity are derived from the upper (lower) bounds on separa- 
tion. 

To develop the definition of p(i,j, k, /), let us define two other separations. 
fl s (i,j) is the single runway separation required for aircraft i following aircraft./'. 
These are the same separations used in the single runway model. fi x (i,k ) is the 
separation required between aircraft i following a departure of aircraft k on the 
other runway. As in the single runway model, these separations are determined 
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from the controller’s point of view, including time to account for uncertainties in 
wind, speed, and position. Let us define as the time of departure of aircraft i. 

Given that aircraft i departs after j on the target runway and k on the other run- 
way, then by definition 

= ma x[tj+n s (i,j), t k + p x (i,k)]. 

In general, the relative values of tj and f* (and hence /,) depend on the unspecified 
history before flight /’ s departure; however, under certain conditions, the separa- 
tion t-tj — i.e., p(i,j,k, l ) — can be computed without knowledge of the prior 
history. 

♦ Markov Property. For any sequence of departures l,j,k such that 
ll s (k,l)< p x (k,j) + p x (j, l ) , 
all prior history is irrelevant in determining 

k k — f^X (^» f) ~b l j 


and 


t,-t j = max[p s ( i,j),p x (i,k) + p x (k,j)\. 

♦ Proof. By definition t. >t l +fi x ( j,l ) , thus 

^ j "b P-x * f) — h "b Px ( ji 0 "b Px (k, y") • 

By hypothesis the right-hand side is greater than p s (k,l) + t,, leading to 
tj + p x (k,j) > (k , /) + tj . 

The two terms above are those whose maximum defines t k , thus the value of t k is 
known in terms of tj. Substituting t } + p x (k,j) for t k in the maximum formula for 
t t , and subtracting tj from all terms leads to the final result. QED. 

Another useful relationship is the following: 

♦ Parallelogram Property. For any departure sequence l,j, k for which the 
Markov Property does not hold, if 

P-s (iJ) + Vx 0\ 0 ^ V s (kJ) + Hx O', k) 
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then 


t, ~tj=Li s (i,j). 

♦ Proof. From the defining maximum formula we note that 
t k -tj = ma x[n s (k,l)-(tj-t,), H x (k,j )] . 

Since tj -t l >fi x ( j,l ) , we have 

t k -tj< max[/x s (k,/) - /J x (j,l), H x (k,j )]; 
and the assumption that the Markov Property is not true leads to 
t k -tj <n s (k,l)-n x (j,l). 

With this result in hand, let us examine the defining relation 
t,-tj = ma x[&(i,j), t k — tj + fi x (i,k)]. 

The second term in the maximum is less than 

(k,l)- fl x F li x (i>k) , 

by the inequality just obtained, and by hypothesis, this bound in turn is less than 
leading to the final result. 

Bounding Separations 

The two properties discussed in the previous section allow direct determination of 
the separation between i and j for some classes k and l. In these cases, we set both 
the upper and lower bound on separation to the known value. For those cases 
where neither property is of assistance, we now describe how to establish bounds 
on the separations. 

The maximum separation between i and j occurs if the prior departure on the tar- 
get runway does not delay flight j by any more than the cross-runway separation 
from flight l. In this case j is leaving as early as possible, considering that flight / 
preceded it on the other runway. If we set tj to the lower bound, t t + H x (j> 0. and 
choose any arbitrary value for t,, then the remaining departure times, including 
can be computed from the defining maximum formulae, and the upper bound on 
the separation between i and j can be computed. 
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The minimum separation between i and j occurs when j is forced to lag l by the 
maximum amount, because of prior history. If 

max_sep(j,l) = max[p s (j ,m) - p x (l,m )] , 


the largest value that tj could take on is t t + ma x_sep(j, l). Assuming an arbitrary 
value for t t and this maximum value tj allows computation of the remaining de- 
parture times and the lower bound on the separation between i and j. 

Both the lower and the upper bounds computed above depend on l and k. The 
bounds independent of l and k are computed by weighted sums of these /, & -de- 
pendent terms. 

Considering More History 

The bounds of the previous section are based on the extreme case for prior history 
These bounds can be refined by explicitly considering prior departure sequences. 
Let us denote the additional flights considered by/,,/ 2 ,/ 3 , ...,/„, each departing 
earlier than the previous one in the sequence. We will use F to denote the entire 
sequence. The flights with an odd index depart from the target runway; those with 
an even index depart from the other runway. The bounds on p p (i,j) are calcu- 
lated as 


bound on p P {i,j) = £(bound due to l,k,F)p xl p xk Y[p f ^ Prf lt • 

l.k.F i=0 z=l 

In practice, we may not need to consider the entire sequence F to bound /— f If 
there is any subsequence f z+2 , f z+1 ,/ z that satisfies the Markov Property, then we 
can determine/., in terms of/ z+] . Given/., and/ z+1 we can determine all subse- 
quent departure times, including the times of interest, /• and t-. Any arbitrary 
value of/ z+1 will do. The capacity algorithm uses recursive code to add history if 
the Markov Property is not true for the last three flights in the current history F . If 
the Markov Property is true, the lower and upper bounds are set to the same 
(computable) value. 

The model user can specify the maximum number of aircraft to add to the history 
F. The larger this maximum, the more accurate the bounds will be, but the longer 
the computations will take. If a particular history sequence has reached its maxi- 
mum size without the Markov Property being true for some subsequence, then 
lower and upper bounds due to the sequence are computed. 
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Before explaining how the bounds are computed, we make the following observa- 
tion: 

♦ Theorem. If the Markov Property does not hold for any subsequence of 
k,j, l, F, then when/„_, is at its earliest time, either all departure times 
within k,j, l, F are based only on same runway separations, or t-tj is in- 
dependent of any further history. 

♦ Proof. Since the Markov Property is not true for any subsequence, the 
cross-runway constraints are not binding on any subsequent flights in 

k,j, l, F when the last two flights in any subsequence occur at their earliest 
times. If additional history requires that some flight /* depart later than its 
unconstrained earliest time, even when/ nl is at its unconstrained earliest 
time — and at this history-constrained earliest possible time for f x , 
t fx + Px (fx - 1 >fx) > */„, + Ps ) , then all departure times after /* 
(including t t , and tf) can be determined in terms of t fs . Furthermore, in this 
situation, adding additional history will not change the relative times of 
departures after f x . If additional history would cause to be later than its 
earlier time, this would cause f x to be deferred by an equal increment, as 
by the assumption it is the accumulated same runway constraints from f n , 
back to f x that have determined t f . A later time for/„_, may also activate 
some other cross-runway constraint, causing /* + , to occur later, but by no 
more than the additional delay to /*; thus,/* would continue to be a point 
from which later departure times can be computed. If there is no such/* 
for the current history, k,j, l, F, this is equivalent to stating that all sepa- 
rations in k,j, l, F are determined by the same runway separations, /u s , 
Q.E.D. 

Now assume that the last flight added is not on the target runway. Then/ n _ ( is on 
the target runway. When /„_, is at its earliest time,y is also at its earliest time. As 
the departure time of /„_, is delayed, it may begin to delay flight j via the accu- 
mulated same runway separations. Thus, the upper bound on separation between i 
and j occurs when t f | -t f + p x (/„_, ,/„ ) the lower bound on t f * ; the lower 
bound on separation occurs when t f * = t fn + ma x_se/?(/ n _ 1 ,/ n ) and the upper 
bound on t, . 

J n- I 

On the other hand, if the last flight added is on the target runway, then/„_, is on 
the other runway. As the departure of f n _ 1 increases from its earliest time, it may 
cause flight k to depart later. The cross-runway constraint between i and k may 
force i to depart later, increasing the time between the departure of i and j. (By the 
theorem, delaying a flight on the other runway either will not change the departure 
time j or will increase the departure time of i and j equally.) Thus the upper 
bound on separation occurs when t f is at its upper bound, and the lower bound 
on separation occurs when t f is at its lower bound. 
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Modeling Current FAA Procedures 

Modeling current procedures requires selecting appropriate values for p s and p x . 
Setting p s (i,j) is described in the single runway model description. We examine 
here appropriate values for p x . 


Departures 


One rule in existing procedures requires a 2-minute departure hold on either run- 
way of a parallel pair separated by 2,500 feet or less after the departure of a heavy 
jet. So if k is a heavy jet, we set p x (i, £) = 2; otherwise, to zero. This rule can be 
turned off by a user input. 

A further restriction when both runways of a pair are used for departures occurs 
when visual separation cannot be applied when a departure is 1 mile from the 
threshold. In this case departures on the parallel runways must be released so as to 
achieve a 1-mile separation. The same departure logic used in the single runway 
model to ensure separation along a single departure path can be used to determine 
the time separation that the controller will apply in this situation. That logic only 
needs to be modified to reflect a 1-mile departure path and 1-mile separation crite- 
rion. 

When ceiling or visibility requires the latter separation criterion to be used, the p x 
value for any pair is the maximum of that required for heavy jet separation and 
that required for the 1-mile separation. 


Arrivals 


Diagonal separation between arrivals to parallel runways may need to be applied 
in IMC. The diagonal separation required depends on the distance between the 
runway centerlines and the radar available to monitor aircraft positions. 

Regardless of the particulars, the diagonal separation can be converted into an 
equivalent separation parallel to the runways, by elementary right-triangle trigo- 
nometry. (The diagonal separation requirement is the hypotenuse; the distance 
between the runway centerlines is one of the shorter sides. The equivalent lateral 
separation is the other shorter side, which can be solved for.) Once the equivalent 
lateral separation is determined, the same procedures used to determine single 
runway controller separations to achieve a miles-in-trail goal can be applied. 
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Appendix F 


Abbreviations 


AROT 

Arrival Runway Occupancy Time 

ASAC 

Aviation Systems Analysis Capability 

ASQP 

Airline Service and Quality Performance 

ATC 

Air Traffic Control 

DTW 

Detroit Metropolitan Wayne County Airport 

FAA 

Federal Aeronautics Administration 

ILS 

Instrument Landing System 

IMC 

Instrument Meterological Conditions 

OAG 

Official Airline Guide 

OASIS 

Information System 

OPSNET 

Information System 

QRS 

Quick Response System 

ROT 

Runway Occupancy Time 

TAF 

Terminal Area Forecast 

VMC 

Visual Meteorological Conditions 
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